Timmerman and Gregorkiewicz Nanoscale Research Letters 201 2, 7:389 
http://www.nanoscalereslett.eom/content/7/1/389 



0 Nanoscale Research Letters 

a SpringerOpen Journal 



NANO EXPRESS Open Access 



Power-dependent spectral shift of 
photoluminescence from ensembles of silicon 
nanocrystals 

Dolf Timmerman* and Tom Gregorkiewicz 



Abstract 

Nanocrystals are widely studied for their tunable optical properties, most importantly increased luminescence 
efficiency and emission energy. Quantum confinement effects are found for many different types of nanocrystals, and 
these introduce a relation between the emission wavelength and the size of nanocrystals. When ensembles of 
nanocrystals with a distribution of sizes are studied, this can have profound effects on their luminescence spectra. 
Here, we show how photoluminescence spectra of ensembles of silicon nanocrystals can shift under different 
excitation conditions, resulting from differences in absorption cross-section of the individual nanocrystal sizes. This 
effect, together with the fact that after a pulsed excitation a silicon nanocrystal can only emit a single photon, 
determines how the distribution of excited nanocrystals changes and leads to the spectral shift for different excitation 
powers. Next to this effect, the influence of different radiative rates in such ensembles is also addressed. These notions 
are important for the interpretation of photoluminescence data for silicon nanocrystals but can be extended to any 
nanoparticle system comprising size-distributed ensembles. 
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Background 

Nanoparticles of many different semiconductor materials 
are intensively studied for their interesting optical proper- 
ties. Quantum confinement in these nanostructures leads 
to strong altering of radiative rates and energy structure. 
In general, stronger confinement in smaller particles gives 
rise to shorter emission wavelengths [1]. Due to limitations 
in size selectivity on this scale, ensembles of nanoparticles 
always show a distribution of sizes and shapes. This results 
in inhomogeneously broadened photoluminescence (PL) 
spectra. For indirect band gap semiconductors, there is a 
special interest in the use of nanocrystals (NCs). In the 
bulk equivalent of these material, electron-hole radiative 
recombination has very low probability, since the differ- 
ence in crystal momentum between the top of the valence 
band and bottom of conduction band can not be compen- 
sated by photon emission. In this case, optical transitions 
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need to be accompanied by phonons, lowering their prob- 
ability. For example, the radiative lifetime of excitons in 
bulk silicon is in the order of milliseconds. In Si NCs, 
the momentum conservation rule is relaxed, as k is not a 
good quantum number any more, and radiative phonon- 
less transitions become more probable. This leads to a 
substantial decrease of the radiative recombination time, 
which can be observed by photoluminescence lifetime 
measurements (from milliseconds to microseconds for 
the smallest Si NCs) [2,3]. 

PL spectroscopy is a powerful tool that is often used to 
determine specific properties of NCs. Most importantly, 
once a relationship has been determined between the size 
and the associated photon energy, it can be used as a rel- 
atively simple experimental way to determine sizes and 
distributions hereof. 

In this work, we study PL spectra of silicon NCs and 
their change for different excitation fluences. A combi- 
nation of size-dependent absorption cross-section and 
saturation effects in these NCs has a profound effect on 
the measured PL spectrum. We show simulations of this 
effect and compare them with experimental data. The 
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effect of excitation conditions on the spectral profile will 
be discussed. Next to this, the influence of different PL 
lifetimes present in an ensemble of NCs is also addressed. 

Methods 

A silicon-rich silicon oxide film with a thickness of 350 
nm was deposited on a quartz substrate by a radio fre- 
quency co-sputtering technique. Subsequent annealing of 
the sample at a temperature of 1,150°C in nitrogen atmo- 
sphere for 30 min induced formation of Si NCs. The layer 
of Si NCs in amorphous SiC>2 created in this way was char- 
acterized by an excess Si amount of 5 vol.% and contained 
NCs with an average diameter of approximately 3.5 nm. 

PL experiments were performed under pulsed excita- 
tion provided by an optical parameter oscillator, pumped 
by the 3rd harmonic of a Nd:YAG laser, with a wavelength 
of 416 nm (3 eV) and pulse width of approximately 9 ns. 
Luminescence was dispersed by a monochromator (Solar 
M266, SOLAR Laser Systems, Minsk, Republic of Belarus) 
and detected by a linear detector (Hamamatsu S 10141- 
1108s, Hamamatsu Photonics K.K., Shizuoka Pref., Japan). 
The entire system was calibrated for its spectral response 
with a calibration source. All measurements were taken at 
room temperature. 

Results 

There is a lot of literature available on the dependence 
of the band gap of the Si NCs on the diameter. For an 
overview hereof, see [4]. For oxidized NCs, it is com- 
monly reported that this dependence follows a D~ b law, 
where 1 < b < 1.5 and D is the NC diameter. The relation 
between emission energy and size of NCs determined by 
high-resolution transmission electron microscopy for the 
sample used in this study has been established earlier [5]. 
It was found that a good fit to the experimental data for 
the optical gap, Eg, is given by the following equation: 

1.86 

Eg(fi) = E Si +-^> (1) 

where E$[ is the band gap of bulk Si, similar as the theo- 
retical dependence determined in [6,7]. When an intrinsic 
yield of PL of 100% is assumed for all sizes [8,9], we 
can obtain a size distribution directly from PL spectrum, 
as the emission wavelength is related to the band gap. 
It is important to note that this is only valid for the 
optically active NCs in the entire distribution. Different 
effects can quench or enhance luminescence from Si NCs, 
which could introduce a size-dependent emission effi- 
ciency. For example, a single dangling bond on the surface 
will quench PL completely [10], and the chance of that 
will be strongly size-dependent. Other effects like exciton 
hopping [11,12] and non-radiative radiation via defects 
[13] can change the wavelength-dependent quantum effi- 
ciencies and, thus, the final emission spectrum. If such 



effects are quantified for a certain system, this could be 
used in order to get information about the entire distribu- 
tion from the determined PL spectra. For the rest of the 
analysis done here, however, we assume that all NCs are 
optically active. 

Rewriting Equation 1 gives the relation between diame- 
ter and wavelength: 

1.86 i 

D(E g ) = (- —)T39. ( 2 ) 

kg — ^Si 

The NCs prepared by the procedure used in this study 
commonly have sizes whose logarithm is normally dis- 
tributed [14-16], described by Equation 3. 

1 (ln(D)-fi) 2 

N(D) = —e ( 3 ) 

Dg \/2tt 

where fi and a are the mean and standard deviation, 
respectively, for the natural logarithm of the diame- 
ter D. By fitting the determined size distribution with 
Equation 3, we can determine these important param- 
eters. Figure 1 shows the result of this analysis, where 
values of \i and a of 1.496 and 0.296, respectively, show a 
good fit with the experimentally determined PL spectrum. 
We note here that the fit actually deviates for larger NCs 
and underestimates its contribution. This is a result of the 
fixed limit that is imposed on the minimum optical gap 
(Equation 1), whereas in realistic systems, this will be less 
strict by participation of phonons in emission. 

In reference [17], it was shown that, under the condi- 
tion that the excitation pulse is much shorter than the 
PL lifetime, every NC can emit a single photon only. 
Multiple excitons localized in the same NC, that have 
been excited by multi-photon absorption, undergo a fast 
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Figure 1 Spectrum and size distribution. PL spectrum obtained 
under high fluence excitation. Inset: The black line indicates the size 
distribution of the ensemble by converting the wavelengths from the 
measured PL spectrum with Equation 1 . The dashed red line shows a 
fit of this distribution with Equation 3 with the fit parameters as 
shown in the figure. 
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non-radiative Auger recombination. Finally, only a single 
exciton remains for radiative recombination. The number 
of excited NCs N* after an excitation pulse can, in this 
case, be described by Equation 4: 

Af*(/ P ump) = N(l - e-Wpump), (4 ) 

where N is the number of NCs in the illuminated volume 
that can emit a photon, a a b s the absorption cross-section, 
and /pump the pump fluence. From geometrical considera- 
tions, it is expected that the a a b s scales with the square of 
the diameter of the NC, which is theoretically supported 
by calculations of optical absorption cross-sections for Si 
NCs of different sizes [18]. We can now write a a b s = yD 2 , 
where y is a proportionality factor. In this case, Equation 4 
can be written as follows: 

Af*(/ P ump) = N(l - e -^W). (5) 

We can use now the expression for the distribution of NC 
sizes from Equation 3 and combine this with Equation 5 to 
arrive to an expression for the number of excited NCs as 
a function of excitation pump fluence and their associated 
band gap energy: 

Jpump) = N(D(E g ))(l - e -^) 2 / P u mp) , (6) 

By converting E g to wavelength with the relation X = jr, 
the associated PL spectrum can be obtained. Qualita- 
tively, this equation describes how the ratio of excited NCs 
for different sizes changes upon illumination with differ- 
ent excitation powers. So, under assumption that the PL 
quantum efficiency is identical for all sizes, it shows how 
the PL spectrum changes with laser power. Larger NCs 
have a larger absorption cross-section, so their emission 
will be saturated at smaller pump fluences than for the 
smaller NCs. This means that, when the pump fluence is 
increased even more, the smaller NCs can still increase 
their contribution to the total (ensemble) emission. As a 
result, the PL spectrum will shift to shorter wavelengths. 
Figure 2 presents PL spectra obtained with Equation 6 for 
six different values of pump fluence of 0.01, 0.03, 0.07, 0.3, 
0.65, and 1. While the intensity increases, saturation sets 
in at longer wavelengths, and the peak shifts to shorter 
wavelengths. A comparison of experimental data and the 
simulations is shown in Figure 3. We note here that car- 
rier multiplication, known to occur in similar systems 
[19], has, in this case, been eliminated by using an exci- 
tation photon energy below the threshold of two times 
the NC band gap. For clarity, the normalized spectra are 
shown. In the simulations, the same ratios of excitation 
fluence are used as in the experiment. The value of y was 
determined from the saturation behavior at a single wave- 
length by the method presented in references [17,20]. In 
this case, Equation 4 is fitted to the power-dependent PL 
intensity to obtain a a b s for a specific wavelength, which 



is subsequently used to deduce a value for y. The sim- 
ulation models the experimental spectra well; the peak 
shifts to shorter wavelengths, and the shift on the short 
wavelength side is more pronounced than that of longer 
wavelengths. The position of the peaks determined from 
the simulation and of the experiments shows a linear 
dependence. The small difference in the absolute numbers 
here is a result of the underestimation of contribution of 
larger NCs and, thus, longer wavelengths, which has been 
discussed earlier. 

Next to a size-dependent absorption cross-section, 
there is also a dependence of the decay time [2] . Typically, 
smaller Si NCs have PL lifetimes in the order of 10 /xs, 
while for large NCs, this increases to a few hundreds /xs. 
This effect leads to red shift of PL spectra for longer detec- 
tion times, which is often observed experimentally [21]. 
For Si NCs, most commonly a stretched exponential decay 
of PL intensity at a single wavelength is observed: 

I(E g , t) PL = I 0 (E g ) • (7) 

where is the stretch parameter; TpL(E g ), the emis- 
sion energy-dependent decay time. In the regime where 
all NCs are excited, we can take the initial distribu- 
tion described by Equation (3) and combine it with 
Equation (7) to arrive to an expression that describes the 
evolution of PL spectrum in time: 

I(E g , t) = N(E g ) • (8) 

For the wavelengths in the range of 600 to 850 nm, we 
have taken PL decays, and the traces have been fitted with 
the stretched exponential function to determine Tpi and 
fi. It was found that ^ = 0.55 over the entire range, and the 
PL lifetime increases exponentially with detection wave- 
length in this region (see inset to Figure 4). For longer 
wavelengths, the xpi was determined by extrapolation in 
order to be able to simulate spectra at different times after 
excitation. Figure 4 shows these results for t = 0, 10, 30, 
50, 80, 100 and 150 /xs. A clear shift of the maximum PL 
intensity is accompanying the change of delay, very similar 
as observed in reference [21]. In the interpretation of such 
a time-dependent behavior, it is, thus, essential to remove 
contributions originating from the above nature. 

Discussion 

The shift of PL spectra to shorter wavelengths when 
increasing the excitation power has been shown to occur 
in different types of NC ensembles. Different explana- 
tions have been suggested in order to account for this 
behavior [22-24], and the importance of the absorption 
cross-section has been noted as well [16,25]. Although 
the above proposed framework does not necessarily fully 
explain the experimentally observed shifts of PL spec- 
tra, it is important to account for the strong influence of 
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Figure 2 Simulated power dependent spectra. Spectra obtained from Equation 6 for pump fluence of 0.01 (red), 0.03 (orange), 0.07 (green), 0.3 
(light blue), 0.65 (blue), and 1 (pink). 



the discussed effects. It should be emphasized that the 
analysis and simulations done here are based on the con- 
dition that the system is excited by a short pulse, with 
respect to the PL lifetime, and that the emerging lumines- 
cence is decaying fully in between the subsequent pulses. 
In experimental practice, this means that the pulse rep- 
etition rate should be small enough to allow all excited 
NCs to return to the ground state. Under continuous exci- 
tation, typical for the majority of experimental reports, 
the PL characteristics change considerably, as equilibrium 
conditions are important in this case. The size-dependent 
decay time could substantially alter the excitation power- 
dependent emission, and an even larger blueshift can be 
expected as the shorter PL lifetime of small NCs allows for 
more radiative recombinations. Now, the question arises 



under which experimental conditions PL spectra should 
be collected in order to get reliable information about the 
physical parameters of the system. For NCs showing long 
decay times, like Si, it is possible to excite the entire dis- 
tribution by a pulsed excitation. After this, all multiple 
carriers simultaneously present in the same NC undergo 
a fast Auger non-radiative recombination, and every NC 
can, in principle, emit a photon. This would then be the 
best way to obtain information about the entire ensem- 
ble. However, high pump fluences might be necessary to 
arrive to this situation, and non-linear and heating effects 
can appear and influence PL properties [26]. Alternatively, 
at low powers where the average amount of absorbed pho- 
tons per NC << 1, the shape of the PL spectrum will 
not change for varying power. This would be the best 
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Figure 3 Experimental and simulated power dependent spectra (normalized). Left: Experimentally obtained normalized PL spectra for 
different excitation powers. The ratio of these powers is the same as the simulations in the left panel. Right: Simulated evolution of the normalized 
PL spectra for ratio's of excitation power identical to that of Figure 2. Inset: Comparison of experimental and simulated peak wavelengths for the 
different excitation powers. 
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Figure 4 Spectral evolution in time. Simulated spectra for different times after a pulsed excitation of the entire ensemble of NCs, acquired from 
Equation 8. The spectral dependence of PL decay time, tpdEg), has been obtained experimentally. Left inset: PL intensity decay for detection 
wavelength of 800 nm.The red line indicates a fit with a stretched exponential function. Right inset: Experimentally determined PL lifetimes for 
different wavelengths (black squares). The red line shows an extrapolation for longer wavelengths. 



region to do comparative experiments, but the low inten- 
sity of PL signal might complicate the practical usability 
of this approach. It is also important to note that, in this 
case, PL will be dominated by contribution from large 
NCs. For more accurate information about smaller parti- 
cles, a scaling based on absorption cross-section should 
be employed. 

Conclusions 

To sum up, it has been shown that PL spectra of typi- 
cal Si NC ensembles show a relatively strong dependence 
on the excitation fluence. A shift of 30 nm in the max- 
imum of PL intensity appears between spectra taken at 
a low power excitation and under saturation conditions. 
It has been demonstrated that this shift can be modeled 
fairly well by taking into account the differences in absorp- 
tion cross-section for different sizes of NCs. This, and 
the notion that each Si NC can only emit a single pho- 
ton after a pulsed excitation, determines the distribution 
of the ensemble that is excited and contributes to PL. The 
evolution of the excited state distribution with respect 
to excitation pump fluence can be modeled fairly well 
by these simulations. Additional simulations, taking into 
account also differences in PL lifetime in the ensemble, 
show spectral narrowing and red shift for peak intensity 
in time-dependent spectra. While the presented modeling 
demonstrates the importance of excitation conditions on 
the PL spectra of ensembles of Si NCs, this approach can 
be extended to all semiconductor NCs, where ensembles 
and size-dependent absorption are imminent. 
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